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Abstract 

The search for the gravitational energy-momentum tensor is often qualified as 
an attempt of looking for "the right answer to the wrong question" . This po- 
sition does not seem convincing to us. We think that we have found the right 
answer to the properly formulated question. We have further developed the 
field theoretical formulation of the general relativity which treats gravity as a 
non-linear tensor field in flat space-time. The Minkowski metric is a reflection 
of experimental facts, not a possible choice of the artificial "prior geometry". 
In this approach, we have arrived at the gravitational energy-momentum ten- 
sor which is: 1) derivable from the Lagrangian in a regular prescribed way, 
2) tensor under arbitrary coordinate transformations, 3) symmetric in its com- 
ponents, 4) conserved due to the equations of motion derived from the same 
Lagrangian, 5) free of the second (highest) derivatives of the field variables, 
and 6) is unique up to trivial modifications not containing the field variables. 
There is nothing else, in addition to these 6 conditions, that one could demand 
from an energy-momentum object, acceptable both on physical and mathe- 
matical grounds. The derived gravitational energy-momentum tensor should 
be useful in practical applications. 
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I. INTRODUCTION 



The notions of energy and momentum play important role in physics M, 0]. These 
quantities are useful because they are conserved. The conservation laws follow from the 
equations of motion, but we can gain important information about the system even without 
explicitely solving its equations of motion. 

For a distributed system (or a field) the densities of energy, momentum, and flux of 
momentum are functions of points labelled by some coordinates x a . These functions combine 
in the energy-momentum tensor T^ u {x a ) 1 that is, the components of T^ u transform according 
to the tensor rule under arbitrary transformations of the coordinates x a (independently of 
whether the space of points x a is endowed with one or another metric tensor). It would be 
embarrassing to use an energy-momentum object which did not transform as a tensor under, 
say, a transition from rectangular to spherical coordinates. Usually, the T^ u is a symmetric 
tensor, T^ v = T u ^. The symmetry of T^ u is required for a proper formulation of the angular 
momentum conservation. The local distributions of T^ v (x a ) are important not only because 
they prescribe some numerical values to the energetic characteristics of the field, but also 
because they can be viewed responsible for the local state of motion of particles and bodies 
interacting with the field. In field theories governed by second-order differential equations, 
one expects the energy-momentum tensor to depend on squares of first-order derivatives of 
the field variables, but not on second derivatives. 

For Lagrangian-based theories, the derivation of the conserved energy-momentum ob- 
ject is closely related to the variational procedure by which the equations of motion are 
being derived (see, for example, [@]). At the beginning it is better to speak about an 
energy-momentum object, rather than a tensor, because at the first steps of derivation the 
transformation properties are either not being discussed or not obvious. In fact, there are 
two routes of derivation. One produces a "canonical" object, and another produces a "met- 
rical" object. The first route takes its origin from Euler and Lagrange. This route does 
not care about transformation properties of the field variables and Lagrangian itself, and 
whether the Lagrangian includes any metric tensor. But what is important is whether the 
Lagrangian contains explicitely (in a manner other than through the field variables) the 
independent variables (coordinates) x a . If such dependence on x a is present, one should 
not expect first integrals of the equations of motion and conserved quantities. If there is 
no such a dependence, some sort of conservation laws is guaranteed as a consequence of the 
equations of motion. 

The second route is associated with the Noether identities. Here one exploits from 
the very beginning the transformation properties of fields and Lagrangians. One requires 
the action to be a quantity independent of any coordinate transformations and, hence, 
one requires the Lagrangian to be a scalar density, that is, a scalar function times the 
square root of the metric determinant. This route produces a "metrical" object, which is 
essentially the variational derivative of the Lagrangian with respect to the metric tensor. 
This object is automatically a symmetric tensor, and it is conserved if the equations of 
motion are satisfied. The conserved tensors are usually understood in the sense that they 
obey differential conservation equations, but one can also derive from them the integral 
conserved quantities if, as is always required, the system is isolated. For radiating systems, 
the fluxes of energy and momentum participate in the balance equations. 
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Both objects, canonical and metrical, are defined up to certain additive terms which do 
not violate equations of motion. These terms are a generalisation of the additive constant 
which arises even in a simplest one-dimensional mechanical problem, when the Lagrangian 
does not depend on time explicitely. It is known that the first integral of the equation of 
motion, which we interpret as energy, can be shifted by a constant. In field theories, the 
additive terms can be used to our advantage. For instance, the canonical object can be 
made symmetric, if it was not such originally, and the metrical object can be made free of 
second derivatives, if it contained them originally. Despite the different routes of derivation, 
the canonical and metrical objects are deeply related. If they are derived from the same 
Lagrangian, explicitely containing metric tensor in addition to field variables, they are equal 
to each other, up to a certain well defined expression calculable from the Lagrangian. 

In traditional field theories, one arrives, after some work, at the energy-momentum object 
which is: 1) derivable from the Lagrangian in a regular prescribed way, 2) a tensor under 
arbitrary coordinate transformations, 3) symmetric in its components, 4) conserved due to 
the equations of motion obtained from the same Lagrangian, 5) free of the second (highest) 
derivatives of the field variables, and 6) is unique up to trivial modifications not containing 
the field variables. There is nothing else, in addition to these 6 conditions, that we could 
demand from an acceptable energy-momentum object, both on physical and mathematical 
grounds. 

When it comes to the gravitational field, as described by the geometrical formulation of 
the general relativity, the things become more complicated. It is often argued that the equiv- 
alence principle forbids gravitational energy-momentum tensor. What is meant in practice 
is that the all first derivatives of any metric tensor g liv {x a ) can be made, by an appropriate 
choice of coordinates x a , equal to zero along the world line of a freely falling observer (along 
a timelike geodesic line). But the first derivatives of g^ v {x a ) can be eliminated along any 
world line, not necessarily of a freely falling observer. And this is true independently of the 
presence and form of coupling of g^ u (x a ) to other fields, and independently of whether the 
9nu(x a ) obeys any equations. Since all components of a tensor can not be eliminated by a 
coordinate transformation, this reference to a physical principle is regarded to be an argu- 
ment against a gravitational energy-momentum tensor, but the argument sounds more like 
a fact from the differential geometry. Despite of this argument, one usually notices that it is 
desirable, nevertheless, to construct at least an "effective" gravitational energy-momentum 
tensor. In practice, this means that we combine some terms of the Einstein equations, in one 
or another manner, into an object, which does not behave as a tensor even under a transition 
from rectangular to spherical coordinates, but which possesses some desirable properties of 
the energy-momentum tensor, and this is why it is an "effective" tensor. And, finally, one 
usually argues that the "effective" tensor becomes the "well-defined" tensor after averaging 
over several wavelengths. Obviously, this transmutation of a pseudotensor into a tensor 
can be done only in an approximate and restricted sense. And, in general, the averaging 
over several wavelengths means that the numerical result will depend on whether we have 
averaged over, say, 3 or 30 wavelengths. 

This shaky situation can be tolerated as long as we are interested only in solving the 
Einstein equations. But this situation becomes risky when we need to know something more. 
It appears that the problem of a rigorously defined energy-momentum tensor may have more 
than a purely academic interest. We have in mind a specific question which was actually 
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one of motivations for our renewed interest to this problem. 

It is likely that the observed || large-angular-scale anisotropies in the microwave back- 
ground radiation are caused by cosmological perturbations of quantum-mechanical origin. 
Cosmological perturbations can be either purely gravitational fields, as in the case of gravita- 
tional waves, or should necessarily involve gravitational component, as in the case of density 
perturbations. To make reliable theoretical predictions one needs to normalize the initial 
quantum fluctuations. In words, this means to assign energy of a half of the quantum to each 
mode. In practice, this implies the availability of a rigorously defined energy-momentum 
tensor for the field in question, which allows to enforce the energy \Two, and not, say, |/uj 
or 30hu, for the initial quantum state. A change in the numerical coefficient would lead to 
the corresponding change in the final results. The preliminary calculations show that the 
contributions of the quantum-mechanically produced gravitational waves and density per- 
turbations should be approximately equal, with some preference to gravitational waves [§J. 
A detailed analysis of the available observational data H seems to favour the gravitational 
wave contribution twice as large as that of density perturbations. Remarkably, the factor of 
2 may turn out to be important when comparing the theoretical predictions with observa- 
tions. This is why, in our opinion, we cannot afford even a numerical coefficient ambiguity 
in such fundamental constructions as gravitational energy-momentum tensor. 

We believe that the difficulty in deriving a proper gravitational energy-momentum tensor 
lies in the way we treat gravity, not in the nature of gravity as such. In the geometrical 
formulation of the general relativity, the components ^(i") play a dual role. From one 
side they are components of the metric tensor, from the other side they are considered 
gravitational field variables. If one insists on the proposition that "gravity is geometry" 
and "geometry is gravity" , then, indeed, it is impossible to derive from the Hilbert-Einstein 
Lagrangian something reasonable, satisfying the 6 conditions listed above. But the geomet- 
rical approach to the general relativity is not the only one available. It is here where it is 
necessary to look at the general relativity from the field-theoretical positions. The general 
relativity can be perfectly well formulated as a strict non-linear field theory in flat space- 
time. This is a different formulation of the theory, not a different theory. The importance 
of looking at theories from different viewpoints was well emphasized by Feynman ||: "if 
the peculiar viewpoint taken is truly experimentally equivalent to the usual in the realm 
of the known there is always a range of applications and problems in this realm for which 
the special viewpoint gives one a special power and clarity of thought, which is valuable in 
itself". 

The field-theoretical formulation of the general relativity treats gravity as a non-linear 
tensor field h tJ,u (x a ) in the Minkowski space-time. In arbitrary curvilinear coordinates, the 
metric tensor of the flat space-time is 7 AU /(a; a ). If necessary, one is free to use the Lorentzian 
coordinates and to transform ^^{x") into the usual constant matrix r}^. The Minkowski 
metric is not an artificially imposed "prior geometry" , but a reflection of experimental facts. 
We know that far away from gravitating bodies, and whenever gravitational field can be 
neglected, the space and time intervals satisfy the requirements of the Minkowski space- 
time. In the presence of the gravitational field, all kinds of "rods" and "clocks" will exhibit 
violations of the Minkowski relationships. This is a result of the universality of the gravi- 
tational interaction (as we understand it today). One is free to interpret the results of the 
measurement as a manifestation of the curvature of the space-time, rather than the action 
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of the universal gravitational field. In this sense, the Minkowski space-time becomes "un- 
observable". But this does not mean that the Minkowski metric is illegitimate or useless. 
On the contrary, it is being routinely used in relativistic astrometry and relativistic celestial 
mechanics. People are well aware of the general relativity and curved space-time. But it 
turns out to be more convenient and informative to store and analyze the data in terms of 
the "unobservable" flat space-time quantities (after subtraction of the theoretically calcu- 
lated general- relativistic corrections), rather than in terms of directly measured "observable" 
quantities. If this is possible and useful in the regime of weak gravitational fields, it can 
be useful for any fields. In fact, for the problem of the gravitational energy-momentum 
tensor, the use of the Minkowski metric allows one to put everything in full order. The 
demonstration of this fact is the main purpose of the paper. 
The structure of the paper is as follows. 

In Sec. II we review definitions of the canonical and metrical energy-momentum ten- 
sors for general field theories. The associated ambiguities, and their relationship with the 
equations of motion, is a considerable technical complication on its own. However, we show 
in detail how the canonical and metrical tensors are related. The main conclusion is that, 
whatever the starting point, the allowed adjustments lead eventually to one and the same 
object satisfying the imposed requirements. We use this general analysis in Sec. IV in 
course of derivation of the gravitational energy-momentum tensor. Sec. Ill is devoted to 
the field-theoretical formulation of the general relativity. We start from the case of pure 
gravity, without matter sources. The gravitational Lagrangian and field equations are given 
explicitely. It is shown that the derived field equations, plus their appropriate interpretation, 
are fully equivalent to the Einstein equations in the geometrical formulation. In Sec. IV, 
being armed with the gravitational Lagrangian and field equations, we apply the general 
definitions of Sec. II for derivation of the gravitational energy-momentum tensor. By differ- 
ent routes we arrive at the energy-momentum tensor satisfying all 6 demands listed in the 
Abstract of the paper. It is shown that this tensor is unique up to trivial modifications which 
do not involve the field variables. We call this object the true energy-momentum tensor. In 
Sec. V we analyze the way in which the true energy-momentum tensor participates in the 
non-linear gravitational field equations. The gravitational energy-momentum tensor is not, 
and should not be, a source in the "right-hand side of Einstein's equations". But it is a 
source for the generalised (non-linear) d'Alembert operator. It is shown that a geometrical 
object most closely related to the derived energy-momentum tensor is the Landau-Lifshitz 
pseudotensor. Their numerical values (but not the transformation properties) are equal at 
least under some conditions. In Sec. VI we include matter fields in our consideration and 
define energy-momentum tensor for the matter fields. The gravitational energy-momentum 
tensor is now modified because of the presence of the matter Lagrangian. However, both, 
gravitational and matter energy-momentum tensors participate in the gravitational field 
equations at the equal footing. Their sum is the total energy-momentum tensor which is 
now the source for the previously mentioned generalised (non-linear) d'Alembert operator. 
The conservation laws for the total energy-momentum tensor are guaranteed by general 
theorems (Sec. II) and are manifestly satisfied as a differential consequence of the field 
equations. The derived equations, plus their appropriate interpretation, are fully equivalent 
to the Einstein's geometrical equations with matter. The final Sec. VII contains conclusions. 
Some technical details are relegated to Appendix A and Appendix B. 
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II. DEFINITIONS OF THE ENERGY-MOMENTUM TENSOR 



Some of the material of this section is known in the literature but we present it in a 
systematic way and in a form appropriate for our further treatment of the general relativity 
as a field theory in flat space-time. 



A. The canonical energy-momentum tensor 

Let us first recall how the notion of energy arises in the simplest case of a 1-dimensional 
mechanical system with the Lagrangian L = L(q, q, t) and the action 

ft-2 



ri2 

S = L(q } q, t) dt . 
Jti 



The equation of motion (the Euler - Lagrange equation) follows from the requirement that 
the action is stationary, 5S = 0, under arbitrary variations of q(t) vanishing at the limits of 
integration (what we will always assume): 

d ±-±( 9 ±\=V. (1) 
dq dt\dq) V J 

The symbol of the total derivative d/dt emphasizes the need to include the partial derivative 
by t if the function L depends on time explicitely. If the Lagrangian does not depend on 
time t explicitely, the equation (|l|) admits the first integral. In this case one has 

dL_dL. + dL.. ^ d_ (dL\ _ d_ (d£\. d_ (9L\.. ^ 
dt dq^ dq^ dt \dq J dq \dq J ^ dq \ dq / 

By multiplying eq. (||) with q and rearranging the terms with the help of (Q), one transforms 
eq. QT]) to 

d ( dL \ 

This equation has the form of a conservation law, and the quantity E = q^r — L, 
called energy, is a constant. With the same success we could call energy the quantity 
E = (q^ — Lj + C, where C is a constant. The equation of motion (|I|) is still satisfied. 

These considerations apply to any field theory described by the Lagrangian L = 
L(qA] qA,a] % a ) where qA(x a ) is a set of variables, and x a is a set of coordinates. The varia- 
tional principle produces the field equations 



dL ( dL 



dq A \dqA, a 



0, (3) 



where the last differentiation with respect to x a includes the partial derivative by x a , and the 
summation over repeated indices is (always) assumed. The field equations are conveniently 
written as ^ = 0, where the variational derivative 5/5 denotes (see, for example, 0) 
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5L(q A ;q Aa ;x a ) _ dL f dL 



5qA dq A \dq A ,c 



(4) 



If Lagrangian depends on second derivatives, the right-hand-side of @) acquires an extra 
term, see Appendix B. 

If the function L does not depend on x a explicitely, one expects that the field equations 
can be transformed into the conservation equations, equal in number to the number of 
coordinates x a . In this case, one has 

dL dL 

By multiplying eq. @ with q Ata , taking summation over A, and making rearrangements 
similar to the ones described above, one obtains, as a consequence of the field equations, 

-6Zl\ =0. 




The expression 

dL 



t£ = q A ,a- €L 



dqA, 



P 



is the canonical (label c) conserved energy-momentum object. The upper or lower positions 
of a and (3 are not essential, but the positions of the first index (a) and the second index 

(/3) are distinguishable. In general, the object t£ is not symmetric in a and (3. 
With the same success we could write for the canonical object 



if the function \l> £ satisfies 



I P - n 4 9L -KPT + ifr P 
oq A .p 



= (5) 



identically or due to the equations of motion (^]). In order to satisfy (^) identically, it is 
sufficient to have ^£ = ip a ^ T T where i/jJ 3t is antisymmetric in (3 and r: = —ipj 13 , so 

that ijj£ T T p = 0. The function ip£ T is usually called a superpotential. By an appropriate 

choice of ^J* one can make the object t£ symmetric in its components. The transformation 

properties of t£ under coordinate transformations are not defined until the transformation 
properties of the field variables and L are defined. 

We now move to covariant relativistic theories. One normally considers physical fields of 
various tensor ranks (scalar, vector, tensor, etc.) in a space-time with some metric tensor. 
The Lagrangian L is required to be a scalar density with respect to arbitrary coordinate 
transformations, that is, L is a scalar function times the square root of the (minus) metric 
determinant. For a better contact with our further study, we consider a symmetric tensor 
field h^ u {x a ) placed in a flat space-time with the metric tensor 7 At[ ,(x a ) written in arbitrary 
curvilinear coordinates x a . The general form for the Lagrangian density is 
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L = LW,h^,h^ ; p) (6) 

where ";" denotes a covariant derivative defined by 7^ and the associated connection 
(Christoffel symbols) C 01 ^. The ^ and C" are functions of x a but they are not dy- 
namical variables, and hence they make the L dependent on x a explicitely. On the general 
grounds, one does not expect the Euler-Lagrange equations to reduce to any conservation 
equations in the usual sense, i.e. in terms of vanishing partial derivatives. However, since 
7 M1/ . a = 0, one can derive a covariant generalisation of the conservation laws, i.e. in terms 
of vanishing covariant derivatives. This is, of course, consistent with our ability to choose 
coordinates x a in such a way that r j^ v will become a constant matrix and C a will all van- 
ish, thus removing the explicit dependence of L on coordinates. Moreover, as we will show 
below, the vanishing covariant divergence will apply to the canonical energy-momentum 
tensor, which is now a manifestly tensorial quantity. 

Let us first give a covariant generalisation to the equations of motion. The action for 
the Lagrangian (|6]) is 

S=-j Ld 4 x , (7) 

where the integral is taken over some 4-volume V. Considering 5h^ v and 5h^ v '. a as indepen- 
dent variations we can write: 

Fit Fit 

6L = -£=-8hr + -^-8hT, T . (8) 

It is easy to check that the operations of variation and covariant differentiation commute. 
Using this property in (Rf) we have 



^-^-jir+i^-shA sir. (9) 



dh^ v \dh^. T / \dh^. T 



Since the quantity Q ^ v is a vector density of weight 1 (i.e. a vector quantity 

times a/— 7) we have 

1 01 5hA = ( J^Sh") . (10) 



\dh> w . T J . T \dh^. T 
Substituting ([5]) into (|7|) and taking into account the equality above one obtains 



5S=- 

c 



d 4 x = . (11) 



At the boundary of integration we have 5h^ u = 0, so the integral of the last term in (|TTD is 
zero. The variations 5h fll/ are arbitrary, and we arrive at the field equations in an explicitely 
covariant form: 

dL ( dL \ _ 



s 



Certainly, one could have obtained the same result in a more familiar way, starting from the 
Lagrangian in the form containing h^ u and the ordinary (rather than covariant) derivatives 
ft/^r (see Appendix A). 

One can now derive the canonical energy-momentum object in exactly the same way as 
was described before. Namely, one multiplies the field equations ([12]) by h^ u - a and rearranges 
the terms to arrive at the covariant conservation law: 

The expression 

^ = tM 7 " t, ^^- 7 "" l ) (13) 

c „ 

is the canonical energy-momentum tensor. We could also define t as 





where \I/ Q/3 is a function such that ^> al3 -p = 0, identically or due to the field equations. The 
conserved canonical energy-momentum object ( pT3|) does not contain second order derivatives 
and is manifestly a tensorial quantity, but, in general, the canonical energy-momentum 
tensor is not symmetric in its components. However, it can be made symmetric by an 
appropriate choice of a non-symmetric ty a/3 . 



B. The metrical energy-momentum tensor 

From the general Lagrangian one can also derive the metrical energy-momentum ten- 
sor. Its derivation relies on the transformation properties of all the participating quantities 
with respect to coordinate transformations. 

An infinitesimal coordinate transformation 

x a = x a - i a {x p ) (14) 

generates the Lie transformations along the vector field £ a , which can be presented as corre- 
sponding variations of the field variables, of the metric tensor, and of the Lagrangian: SW, 
5 r y liu , and 6L, respectively. Since the Lagrangian @ is a scalar density, its variation is a 
total derivative 

5L = (LC), a ■ (15) 

The change in the metric tensor is 

fly*" = _f w _ f . (16) 

And there is also a corresponding change in the field variables 
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5h aP = ^h a/3 . a - h aa f. a - h^i a . a (17) 

but we will not need to know its concrete form for this derivation. 

Taking into account (|1|) and assuming that the vector field £ a vanishes at the boundary 
of integration, we conclude that the variation of the action must be equal to zero: 

SS = ~ c J 5Ld' l x = . (18) 

On the other hand, we know that an arbitrary variation of L, not necessarily caused by ([TiD, 
has the general form 

5L = -^-5hT + -^-S-f + A a , a (19) 

where 

a* = -^-Shr + -j^-tr ■ 

In writing this formula we took into account the fact that the first derivatives of 7 Miy partic- 
ipate, through the Christoffel symbols, in the covariant derivatives of the field variables. 
At this point we require that the field equations are satisfied 

51 (20) 



8h^ 



so the first term in flT9"|) is zero. In the second term of flT9"|) we use the specific variation 
Then, eq. ([18]) acquires the form 

(2i) 

where we have also used the following equality: 

9L dL 

-ImlvpT, • (22) 



The second integral in (|2~ID transforms into a surface integral and vanishes under appropriate 
boundary conditions for £ a . Since the functions £, p {x a ) are arbitrary, we finally obtain: 



5L 



. (23) 



Pa , 



The metrical (symbol m) energy-momentum tensor t ^ is defined as 

m in, 2 5L . , . 

t 1 ™ = --=-. , 24 

so that eq. (^) takes the form of the covariant conservation law (valid only on solutions to 
the equations of motion (p0|) ): 
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(25) 



As before, one can also write for the metrical energy-momentum tensor 

T^ = 1=-^ + $^, (26) 



where the function ^ v satisfies <& Mt = identically or due to the field equations. The 
derived conserved object (|24]) is automatically symmetric and a tensor, but, as a rule, it 
contains second order derivatives of the field variables, even if the Lagrangian does not 
contain them. They are generated by one extra differentiation in the definition of the 
variational derivative (see the second term in eq. (|])). However, by an appropriate use of 
and the field equations, all second derivatives can be removed, as we will discuss in 
detail later on. 



It is important to note that nothing in the derivation of eqs. (p4[), (|2q ) actually required 
the 7^ to be a metric tensor of a flat space-time, that is, to have the curvature tensor 
constructed from 7^ equal to zero. One can still formally arrive at the equation, similar 
in structure to eq. (p5j ), in arbitrary curved "background" space-time, where the covariant 
derivatives are now being taken with respect to the curved metric. This is an element of a 
field theory in the "background" space-time, which is useful in some applications (see, for 
example, M, ||). As soon as the field equations are satisfied, the corresponding covariant 
"conservation laws" must be valid. However, in this case, there is not and there should not 
be, in general, any conservation laws in the usual sense. First, one normally encounters 
severe integrability conditions for the field equations. The number of independent solutions, 
in the sense of the Cauchy problem, can be diminished, or the solutions may not exist at all. 
Second, the vanishing covariant divergence cannot be converted into the vanishing ordinary 
divergence. This is a well known formal obstacle, but it has deep and clear physical reasons: 
the "background" space-time is by itself a gravitational field which interacts with a system 
and can exchange energy with the system. For instance, even in the simplest Friedman- 
Robert son- Walker space-times, gravitational waves can be amplified and gravitons can be 
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Returning to the strictly defined energy-momentum tensors, we will now show that the 
canonical tensor fllTf) and the metrical tensor (E3) are closely related. 



C. Connection between metrical and canonical tensors 

The metrical tensor (|24D and the canonical tensor (|13D are derived from the same La- 

,1— k m c 



grangian (P), so one expects them to be related. To find the link between t ^ and we 



m 



return to the derivation of t ^ v based on the infinitesimal transformation fll4]) . 
It is convenient to write the variation ([15]) in the form 

SL = {LC\a ■ (27) 

The replacement of the ordinary divergence by the covariant one is allowed, because the 
differentiated quantity {L^ a ) is a vector density. We can also write the general variation 



(|T9|) in the form 
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5L 



5L 



Sy^ \dh a P. T , 



+ 



' dL 



$1 



(28) 



In writing this expression we took into account ([LC]) and the fact that 

dL 8L 

We will now show that the differentiated quantity in the last term of fl28"|) is also a vector 
density, so that the ordinary divergence can be replaced by the covariant one. Indeed, from 
the structure of (|5]) it follows that 



dL 



dL dh^, p _ 2 dL h „ p dC% 



Since 

dC 



dy 



a/3 



97"" T dh a P. tP dji u/ dh^.p dy 



-aW-TPX + ^^7oA + WIIPp + ^y ap ~ 7™7p«7/ 3 a - 7™7*07«a) (29) 



it is now clear that the quantity 9 ^ 5y^ u is a vector density Thus, we can rewrite (p8|) as 



SL = ^SV* + + ( 2 



9L h *p dC % 5 , 9L ^ 
dh<#. p dr v T dh^. . 



(30) 



The expression (p0|) is valid for arbitrary variations, and hence it is valid for specific 
variations ([RJ), fllTD caused by (fH]). Therefore, the difference between (|27| ) and ( j30|) must 
be equal to zero. Substituting ( |TTD and (|TBp into this difference, and combining in separate 
groups the terms which contain ^ <J . r and £ a . T .\ one obtains the equality which should be 
true for arbitrary vector field £ a (x a ) : 



5L 

8h a P 



' a+ \dh^, T 



h ap . n - 5ZL 



-2-^V + 



dL 



5L 



dh a P. T 



dh a P.p dy au x 



dL 



dC\ 



dh^.p dy au 



-7 



dL 



dh aa . 



(31) 



The coefficient in front of £ CT is identically zero, because all the terms cancel out. To 
check this one has to recall the definition of the variational derivative 



5L dL 



dL 



Sh a[) dh a(3 I Q h af3 
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to use (|A4|) and ( |X3| ) for -J^$ and ) > an d to take into account 



dL h a P 4- ^ h a P - T 

a ;<r T OI — '1 ; T ;o- — 



The last term in (0), which contains £ CT .;v T , is also identically zero. This is true because 
the ^' cr . A . T . is symmetric in the indices A, r whereas the coefficient is antisymmetric in these 
indices. To show this in detail, we denote this coefficient ^— r )ipJ x and rewrite it using 
formula (^D|): 



dL 



h* p 



dC\ 



dL 



dh&.p dj 



-7 



dh c 



-h 7 



dL 



dh°P., 



-h x ?- 



dL 



■K* + 



dL 



dh a ?., 



14» 



dL 



dh a P. 



-T' - 



dL 



-7 



aX 



(32) 



It is now clear that ipj x = —ip a XT ■ So, we are left only with the term which contains £ CT . r . 
Since the vector field is arbitrary, this gives us the equation 



-2- 



5L 



$7 



pa 



jPr + h a/3 



dL 
dh a P. T 



-51L \ - 2 



5L 



h™ + V=7^ A ;A = • 



5fc 



Using in the first two terms the definitions of t ^ v and t^ v and formula (|22|), we arrive at 
the universal relationship 



fiVT 



1 -_Y a w^ = o. 



/_ 7 



(33) 



Assuming that the field equations are satisfied (the last term vanishes) we can finally con- 
clude that 



t 



Thus, the metrical and canonical tensors are related by a superpotential whose explicit 
form is given by eq. ([32]). (This derivation is similar to the one given in |JL1[|.) Obviously, 
the conservation laws are satisfied because ^ VT . T . v = . 



III. FIELD THEORETICAL FORMULATION OF THE GENERAL RELATIVITY 

The field theoretical approach to the general relativity treats gravity as a symmetric 
tensor field h^ v in Minkowski space-time. This approach has a long and fruitful history. In 
fact, in the early days of special relativity, Poincare and Einstein himself started from an 
attempt to give a relativistic generalisation of the Newton law. Even after the acceptance 
of the geometrical viewpoint, various aspects of this approach have been worked out in 
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numerous publications [0- |[28|| , 0, |19| - to name only a few. (One can also find references 
20fl useful.) We will follow a specific scheme developed in || and ]l{J, as a continuation of 



the line of reference p8 



The gravitational field h^ v (x a ), as well as all matter fields, are defined in the Minkowski 
space-time with the metric tensor 7^^(x Q ): da 2 = r ) iliV dx^dx v ' . The matrix 7 Ml/ is the inverse 
matrix to -y^ u , that is, 7 o/3 7/3i, = 5™, and 7 is the determinant of the matrix 7^. The 
raising and lowering of indices are being performed (unless something different is explicitely 
stated) with the help of the metric tensor 7^. The Christoffel symbols associated with 7^ 
are denoted by C T ^ V) and the covariant derivatives are denoted by a semicolon The 
curvature tensor of the Minkowski space-time is identically zero: R a p^u{l pa ) = 0. 

In terms of classical mechanics, the field variables h^ v are the generalised coordinates. 
Their derivatives W-r (a third-rank tensor) are the generalised velocities. It is also conve- 
nient (even if not necessary) to use the generalised momenta P a canonically conjugated 
to the generalised coordinates h^ v '. The object P a is a third-rank tensor, symmetric in its 
last indices. We will also need a contracted object P a = P T aT = 5 v P^ ail . 

The use of W and P T ^ U as independent variables is an element of the Hamiltonian 
formalism, which is also known as the first order variational formalism. We will start from 
this presentation, and then will consider the presentation in terms of and W.^. It will 
be shown that the derived field equations are fully equivalent to the Einstein equations in 
the geometrical formulation of the general relativity. 



A. Gravitational field equations in terms of generalised coordinates and momenta 



The total action S of the theory consists of the gravitational part S 9 and the matter 
part S m : S = S 9 + S m . We will include the matter part in our consideration later on (Sec. 
VI). The action for the gravitational field is 



S 9 



d 4 x 



where the Lagrangian density L 9 is 



L 9 



"7 



2k 



(34) 



and k = 87rG/c 4 . It is now clear that the quantities P T ^ U are indeed the generalised momenta 
because 



-7 



P T 



dL 9 



dL 9 



2k " ^ dh^. T dh^\ T 
The tensor P£ v is related with the tensor K T ^, U originally used in || by 



p\ v = -k\ v + -s;k u + -6;k, . 

To make the part of L 9 , which is quadratic in the momenta P^ u , more compact, we will 
also write the Lagrangian in the equivalent form: 



14 



V 



-1 



2k 



i^pa pa _ pr nu 

\<* r per 2 WT P a a f- 



(35) 



where 



(36) 



and n^ a(S WT = n v ^ WT = = fi a/3/ ^™ • 

The gravitational field equations are derived by applying the variational principle to 
(|34|) and considering the variables h^ u and P" as independent. In this framework, the field 
equations are 



dU ( dU 



and 



dL9 



dL 9 



dP> 



dP a 



(37) 



LIU- 



Obviously, the term d p]f in fl3~T|) is zero for the Lagrangian fl34"|). Calculating the deriva- 
tives directly from (|34]) and introducing the short-hand notations for the corresponding 
expressions, one obtains 



2k 5L 9 



— r =_P a _P Q P/ 3 _i pp_n 

' (UV 1 pz,;a J ^/J- 1 j/ Q T g -1 ^t-l jy <J 



(38) 



2k 



fj.ii 



(39) 

Using the f2-matrix introduced above we can rewrite eq. ( |39"D in the compact form: 

hr. T = ^ a ^ T p% . (40) 

Equations r Mi , = and f T ^ u = form a complete set of equations in the framework of the 
first order variational formalism. 



B. Field equations in terms of generalised coordinates and velocities 

We will need the field equations in terms of the gravitational field variables h^ u and 
their derivatives. We will derive the equations from the Lagrangian (|35| ) written in the form 
containing the generalised coordinates and velocities. This is an element of the Lagrangian 
formalism, known also as the second order variational formalism. To implement this program 
one has to consider P a as known functions of h^ v and h^ u . a and to use them in the 
Lagrangian (p5|). 
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The link between h) iV and P T ^ U is provided by eq. (f40l). To solve equations (^) with 

respect to P T „ V we introduce the matrix fiZjL, , which is the inverse matrix to Vt a(3 ^ v TU and 
satisfies the equation 



r\pva/3 q — 1 1 " l l J 
14 UT^ L prJ pv 



(41) 



The explicit form of the Q 1 -matrix is not needed for the time being, but it will be given 



below where required. We will only use the symmetry properties of the i\ J pcr ™ which are the 

-l TU = Q-l TUJ = Q-l TU = Q- 

"pvpa vppa fiver p papv 



same as the symmetry properties of the f2-matrix: Q 



_1 UJT 



By multiplying the both sides of eq. (pLOf) with 



i-i T V> 



papv 



one obtains 



h pcJ 

* L pcjpv "' ; 



Now we substitute eq. 



pv papv ,ul ■ 

(42) 

into eq. fl35|). The Lagrangian takes the elegant form 



LP 



4k s W ft ^ n > 



(43) 



which is manifestly quadratic in the generalised velocities h pv ;r . The dependence on the 
generalised coordinates (as well as on the metric tensor 7^) is contained in the Q^ 1 
tensor. The Lagrangian L 9 belongs to the class of Lagrangians (§) studied in Sec. II. 
The field equations in the framework of the second order variational formalism are 



5L9 



-(- 



dL 9 



Sh^ dh» v \dh^ u - T 



. 



In more detail, we have 



QQ-l " T 

U ^paal3 h pg , a f3 _ r, ( Q-l WT, a/9 \ 



. 



The first term can be calculated by differentiating (|4l|) with respect to h pu and taking into 
account (|3~6"D. This gives 



«q-1 ^ 

Q h pv 



StSt - -St St 



Q-l T7r Q-l W£ 1 Q-l T7r Q-l We 



The second term requires to recall the rules of the covariant differentiation applied to the 
tensor which, in turn, is a function of h pu and r ) pu : 



dhP a 



Combining all together, one arrives at the field equations which are manifestly the second- 
order differential equations in terms of h pu : 



n 



-1 

pva.fi 



h af3 ._ ur _ T - stst - ^stsi 



9 Q-1 T7r Q-l We _ Q-l T7r Q~l U)t 
ZlL ptppv lL a<j>af3 lL ptppa lL v<j>af3 



h pa . T h al3 ]UJ = . (44) 



Certainly, one arrives at exactly the same equations by substituting P T found from eq. 



(see eq. fl42|) ) directly into fl38|) 
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C. Equivalence of the field theoretical and geometrical formulations of the general 

relativity 

We will now show that the entire mathematical content of the general relativity (without 
matter sources, so far) is covered by the Lagrangian fl34|) , or by its equivalent form (fj3[). 
We will demonstrate the equivalence directly at the level of the field equations, rather than 
at the level of the Lagrangian ( |34|) and its Hilbert-Einstein counterpart. The derived field 
equations (j38|), (|39| ) can be rearranged by identical transformations into the usual Einstein 
equations. 

First, we introduce a new tensor field g^ v (x a ) according to the rule: 

V=gg» v = V=i(Y u + JO (45) 

where g = det\gn U \ and the tensor g^ u is the inverse matrix to the g^ u matrix: 

Q^Qua = K ■ (46) 

Let us emphasize again that the tensor g pv is the inverse matrix to g^, and not the tensor 
g^ v with the lowered indices, g^ v ^ Ifialvpg^ ■ For the time being, we do not assign any 
physical interpretation to the tensor field g^, we only say that the functions g^ u {x a ) and 
gixu(x a ) are calculable from the functions h^ u {x a ) and 7 M!/ (a; a ) according to the given rules 

(H), (11). 

The introduced quantities allow us to write the f2-matrix as 



QpO-Qf/3 _ V 9 

J L UJT „ , 



We can also give the explicit form for the ^KlpJ ■ By multiplying the both sides of (f4T|) 
with l^[25i(5%5l + 5%5l)g e/3 - g 5u {2g a4> g px - gapg^x)] one obtains the explicit form of the 
f^-matrix: 

%l p r = + WIW P % + Kfyg.x - g^{g w g v . + g vp g^ - g^g p(T )\ . 

We now want to calculate the quantity T T defined by the expression 



I 

T 



r V = ng rX (g^M + gxw - g^,x) • (47) 



By replacing the partial derivatives with the covariant ones we get 

rv = c \u + \g TX {g^u + gx^ - g^x) ■ (48) 

Now we want to trade g pU]T for [ y /—gg^ ll '). T in order to have quantities easily expressible in 
terms of 7^ and h a/3 . By differentiating (|^) one obtains 
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9pv; 



(49) 



Using the formula for the differentiation of determinants, we can write 



-99 



pa\ 



:9af39 



pa i 



'99^ 



Substituting (||) and Q50|) in (|48|) we obtain 



(50) 



r r = c r a 

pv ^ pv > 2 



-g 



-K9pp{^99 Pa \v - ^9 V p{^T99 pa \p + 9 TX 9pp9uAV^99 Pa );X+ 



1 



"7 



pupa 1 g M TV pa 1 g i/ Tppa 



-99 pa )-,x 



Finally, taking into account {^/—gg 



pa\ 



'—^h pu . a and recalling (^), we arrive at 



-pr _ sit _ pr | ^rr p i p 
L pv ^ pv 1 pv ~ ^p 1 v 2 u v 1 P • 



(51) 



Now we want to use (plQ and calculate the quantity _R Mi , defined by the expression 



p _ pa 

1 /if , 



^■pa ^pa _j_ pa p/3 pa p/3 

2 pa,v 2 wo,/i /if a/3 A*/3 



(52) 



The C T ^ part of r r Ml/ produces a series of terms which combine in the Ricci tensor R pi/ of 
the flat space-time. The ordinary derivative of the tensor PL, plus all the terms containing 
the product of P% with C^ av combine in the covariant derivative of PL,. All other terms 
produce quadratic combinations of Pfy. In the result, we arrive at 



R 



pv ± L pv 



P — ( P a i pa p P p p 

lh uv \ 1 fiv;a * 1 p/3 1 va ^ 1 u 1 i 



p ± v 



(53) 



Since = we conclude that the field equations (|3lf) are fully equivalent to the equations 



Rpv — . 



(54) 



The remaining step is the matter of interpretation. We can now interpret the quantities 
g a p as the metric tensor of the curved space-time: 

ds 2 = g^dx^dx" . (55) 

Then, the quantities (f|7]) are the Christoffel symbols associated with this metric, and the 
quantities (|52"D are the Ricci tensor of the curved space-time. Finally, equations fl54|) are the 
Einstein equations (without matter sources). 
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IV. THE GRAVITATIONAL ENERGY-MOMENTUM TENSOR 



Being armed with the definitions of the energy- momentum tensor (Section II), as well as 
with the gravitational Lagrangian and field equations (Section III), we are now in the position 
to derive the gravitational energy-momentum tensor. We will derive both tensors, metrical 
and canonical, and following the general theory of their connection, we will find explicitely 
the superpotential which relates them. We will show that the requirement that the metrical 
tensor does not contain second derivatives, and the requirement that the canonical tensor 
is symmetric, produce one and the same object which we call the true energy-momentum 
tensor. This object satisfies all the 6 demands listed in the Abstract of the paper. 



A. The metrical tensor 



The metrical energy-momentum tensor defined by eq. (|23 ) and derived from the La- 
grangian density fl34[) has the following form: 



ml 1 1 

Kt r = - 7 ^%P% + [VV - fTW + hn\{P a pP P L ~ + V , (56) 



where 



Q 



Expression (^) was obtained by direct calculation of the variational derivative (or alterna- 
tively, see Appendix B) 



5L 



dL 



dL 



\tt> 



Pi, 



(57) 



and no further rearrangements have been done. Obviously, tensor fl56[) is symmetric in 
its components, but it contains second order derivatives of h pu which enter the expression 
through the Q pv term. We want to single out the second derivatives of h pu explicitely. 

By making identical transformations of the Q pu term one can show that the Q pu contains 
a term proportional to r M „ and terms proportional to f T pv and its derivatives. All these terms 
are equal to zero according to the field equations ( |58"D and Q59]). After removing these terms, 
the remaining expression for Q pu is as follows 



Q 



1 



(TC + W) 



2 V p 



af3pa(piT p\ _p p \_ a.Ppa pr 



(58) 



where 



a(3pa 



h *ayf3 + frpayf} + ffdrf* + h pP Y a + h aa h pl3 + h pa h al3 - h pa {^ a(3 + h af3 )] . (59) 



The remaining expression (|58|), together with other terms in (|56|) , reduce the Kt pv to 
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m 



r a/3 



1 



hh va , a h^. iP - hf a . p h vp . a ) + i(-2/« a/3 + bTh^ + h va h^), a , p (60) 

where the subscript \ r indicates that the energy-momentum tensor was reduced on the 
equations of motion. 

The last group of terms in fl6"0|) still contains second order derivatives of ft/"', but they all 
can be removed by a special choice of superpotential. Indeed, the symmetric function <&^ u 
participating in ( P6"D and satisfying ^ v . u = can be written as 

= (0^ + <f)^ al3 ). a .p (61) 

where 

To remove all the second order derivatives, we require 

^{~2h^h af3 + h^h Ul3 + h Va h^), a] p + (0^ + 0^ a/3 );a;/3 = 0. (63) 

The unique solution to this equation (up to trivial additive terms which can possibly contain 
7 M1/ but not the field variables ft/ 1 ") is: 

= _^h^ - h av h^). (64) 

With the help of the superpotential (§4]), we can now cancel out the terms j(—2h fiu h al3 + 
h m h uf3 + h ua h^) W f3. The remaining part of (^) does not contain any second order deriva- 
tives at all. To write the remaining part in a more compact form, we replace the generalised 
momenta by the generalised velocities with the help of (0), and use the shorter expressions 
g a/ 3 and g a/3 according to their definitions (^) and (pEB]). As a result, the metrical energy- 
momentum tensor (0), transformed with the help of the field equations and an allowed 
superpotential, takes the following explicit form: 

KIT = -^W v , p W° , a - 2hT, a h^. 3 + 2g^g a ^, (J hr, p + g^g ap h^, a h^, p - 
W a gp P h vP ;«h p °., a - 2g va gfip h^, a h^, a + 

i(2<?'V - g^g^g^p - g^g^,^,^ (65) 

where g a p and g af3 are short-hand notations for the quantities ([45|) , ([4~6|). This object is a 
tensor with respect to arbitrary coordinate transformations, symmetric in its components, 
conserved due to the field equations, free of second derivatives of ft M ", and unique up to 
additive terms not containing ft/ 1 ". This derivation required the use of an allowed superpo- 
tential. The last step is to show that the energy-momentum tensor (|65D can also be derived 
according to the original definition p4"|), without resorting to the use of a superpotential. 
The tensor (|65D will be derived from a modified Lagrangian, which produces exactly the 
same field equations as (|38D and (^). This is what we will do now. 
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B. The constrained variational principle 



Let us write the modified Lagrangian in the form 

(66) 



2k 



where Rappa is the curvature tensor constructed from 7^. Obviously, we have added zero to 
the original Lagrangian, but this is a typical way of incorporating a constraint (in our case, 
Rap/3a — 0) by means of the undetermined Lagrange multipliers. The infinitesimal variation 



(HID of the metric tensor (and even its exponentiated finite version) do not change the 
condition Rapp^ = 0. The multipliers A a/3pa form a tensor which depends on and W 
and satisfy 

J^aPpa _ _J^pPaa _ _^aupfi _ ^Pacrp 

The variational derivative of A Q|8pa R ap pa with respect to the metric tensor 7^ is not zero, 
and therefore the added term will affect the metrical energy-momentum tensor. However, 
the added term does not change the field equations, since the variational derivative of this 
term with respect to the field variables h pv will be multiplied by the R ap pa and hence will 
vanish due to the constraint. 

The metrical energy-momentum tensor fl24]) directly derived from ( |6~6"1) is now modified 
as compared with fl56"l): 



ml 1 1 

Kt He = -Y u h^, a p a pa + [Y P Y° - ^ u (Y a + hn}(P a P pPL - ^PpP.) + 

_ ^af3 + A^ a ?). a .p (68) 

where the subscript | c indicates that the Lagrangian (|66|) has been used. The entire modifi- 
cation amounts to the last two terms (with double derivatives) in (|5"8"D, which immediately 
suggests its connection to modifications at the expense of superpotentials flST]), (|62|). (For 

m 

a detailed derivation of the last two terms see Appendix B.) As before, the tensor Kt^ u \ c 
contains second derivatives of h pv in the Q pv term. But the originally undetermined multi- 
pliers k a(ipa will now be determined. They can be chosen in such a way that the remaining 
second derivatives of h) w (which could not be excluded at the field equations) can now be 



removed. The equations to be solved are similar to equations (|63|). Their unique solution is 



A puaP = _h h *P h pv _ ffwtflt). 
4 

Thus, the energy-momentum tensor (|65|) satisfies the last remaining demand: it can be 
derived in a regular prescribed way (El) from the Lagrangian (|66|). 



C. The canonical tensor 

The gravitational energy-momentum tensor (53) satisfying all the necessary demands 



has been derived along the "metrical route". We will now show that the symmetrisation 
procedure of the canonical tensor leads to the same object (|65|). 



21 



directly calculated from the Lagrangian 



The canonical energy-momentum tensor 
density fl43[) has the form 

1 

4 1 

It is convenient to use here and below the quantity P T ^ U as a short-hand notation for 
^^laf^ T ^ 3 -,uj m agreement with fl42|) . Then, the re t"" takes the compact form 



Kt 



pi' 



(69) 



As expected, the canonical tensor is not symmetric. It can be made symmetric (see II A) 
by an appropriate choice of ^/ fJ-u . We will do this on the basis of the universal relationship 
(|33|) between the (symmetric) t^ u and the (non-symmetric) t^ v . 
The relationship in question is 



- ret ^ + ret"" + + 7 MQ ^r a/ 3 = 

where ip^ UT is calculated from the Lagrangian (33) according to (P2f): 



(70) 



Kip^ LUT = -[P^{Y a h u/3 - -fh 1 *) + P^ aT h^ - 7 a ^ T/3 ) + Pap(l aT h^ - -ft??)]. (71) 

As any superpotential does, the tensor tp tJ-UT satisfies the requirement ?/> A " yT . T . u = 0. One can 
easily check the validity of (|70"D if one combines (|56"1), (|69|), (f7T|), (|38|) and uses the following 
identities: 



1 



(7^ + ^)(V; 



1 



P a Pp ) = h V(T . n P p , 



;p eta 



pa i 



2 i" pa 



We now assume that the field equations are satisfied. The last term in (|70f) drops out. 

The metrical tensor t ^ v reduced at the equations of motion is given by (|60|) . We need 
also to reduce the third term in ( ffO| ) at the equations of motion. First, we differentiate the 
expression (|7lD and make identical transformations to rearrange the ip flUT . T : 



Kip 



pUT 



h^(Y a + h ua )- -/^( 7 Q/3 + h a(3 ) 



■(-2h^ u h a/3 + h^h" + h ua h^), a . p 



(72) 



where 



X 



pv 



h^(Y a + h ua ) - -h^^ + h al3 ) 



h^(Y a + h va ) - -h^^ + h a/3 ) ( 

2 . ;r \ 



pp pa 

r aa r I3p 



^{h va , a h^, p -h' ux ,ph vP , a ) . 
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Now we drop the term proportional to the field equations. The remaining part of ipi^ 7 ' . T is 

K^ UT . T \ r = + ^{-2W v h aP + W a h v(i + h va h^). a . p . (73) 
On the field equations, the original relationship (|70|) reduces to 

7^| r =t^+ V VT , T \r- (74) 

The last term in expression ( [73] ) for . T \ r cancels out with exactly the same term in 

expression (|60|) for t^ v \ r - After this cancellation, what is left on the left hand side of ( |74|) 
is the metrical tensor t^ u described by formula (po]). On the right hand side of ([T3]) we will 

c c 

get a symmetrised (subscript \ s ) canonical tensor Kt^ v \ s = Kt^ v + x^ v ■ Note, that since 
ip^ . T . u = 0, it follows from (|72|) that x^-u = on the field equations r af3 = 0. Thus, we 
arrive at the equality 

Since the metrical tensor ( |65"D satisfies all the 6 demands listed above, and since it can also 
be obtained as a result of symmetrisation of the canonical tensor, we call it true energy- 
momentum tensor (and write it without any labels or subscripts). 



V. GRAVITATIONAL FIELD EQUATIONS WITH GRAVITATIONAL 
ENERGY-MOMENTUM TENSOR 



We have derived the gravitational (true) energy-momentum tensor (|65|) from the gravita- 
tional Lagrangian L 9 according to the general definition (0). We know that the conservation 
laws t^ v . v = are guaranteed on solutions to the field equations. The non-linear nature of 
the gravitational field W makes the field a source for itself. The question arises as for how 
the Kt^ u participates in the field equations. To answer this question one needs to rearrange 
the field equations and single out the nt^ v explicitely. One can proceed either from equations 
(|38|) or from equations (|44|). A simpler way is to take the following linear combination of 
the field equations (|38|): 



a/3 







(75) 



and to use the link (|4^) in order to exclude P a ^ v - After putting all the terms in a necessary 
order, the field equations (|75|) take the following form 

I[(y- + + h af> ) - {Y a + /Ofr"" + h v % a .p = - 



2h^ a , a h^ ]f3 + 2g^g^h v(} , c h^, p + g "* ' g ap h^ ;a h»° ;/3 - 2g^g Pp h^ ;a h^. a 



2g™gp p h^. (T ho° ;a + -(2g» 5 g™ - g^g^){2g pa g^ - g aP g p<7 )h^. s h^. >u 



(76) 
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On the right-hand side of (|76|) we have exactly the energy- momentum tensor (|65| ), so the 
field equations can be written 

+ W v ){ 1 a P + h aP ) - (y° + h>"*)(Y P + h vP )}, a , p = KIT. (77) 

The left-hand side of this equation is the generalised differential wave (d'Alembert) operator. 
So, the gravitational energy- momentum tensor Kt^ v is not, and should not be, a source term 
in the "right-hand side of the Einstein equations" , but it is a source term for the generalised 
wave operator. 

Replacing the sum (7^ + h^ v ) by the shorter expression g^ v according to the definition 
([45|), we can also write the gravitational field equations fl77l) in the form 



\[(-9)(gr<f* - g m g v %^ = kT . (78) 



We know that the field equations (|7q) are fully equivalent to the Einstein equations (see 
Sec III). In the geometrical approach to the general relativity one interprets the quantities 
g^y as the metric tensor of a curved space-time (|55|). It is interesting to ask if there exists 
an object in the geometrical approach, which would be somehow related to the energy- 
momentum tensor t^ v derived here. [The description of energy in the general relativity is, 
of course, a matter of a long time effort by many people who used different approaches. We 
would like to mention at least some of works |21||- ||32|| , which influenced our understanding 



of the problem.] For this purpose we use for the first time the available coordinate freedom 
and introduce the Lorentzian coordinates. This means that the metric tensor / y fJiU (x a ) is 
being transformed by a coordinate transformation to the usual constant matrix 77^. In 
other words, one makes 700 = 1, 7n = 722 = 733 = — 1, the rest of components zeros, and the 
determinant 7 = — 1. Then, all the covariant derivatives can be replaced by the ordinary 
ones, and all the derivatives of the metric tensor 7 Miy vanish. Writing the expression flSo] ) 
for in Lorentzian coordinates (subscript \l) and using quantities g^ v instead of bT one 
finds that 

n L = (-g)ttl 

where t^ L is the Landau-Lifshitz pseudotensor 0. The field equations ( [78] ) written in 
Lorentzian coordinates take the form 



So, the object most closely related to the derived energy-momentum tensor t^ u is the Landau- 
Lifshitz pseudotensor t^ L times (— g). Their numerical values (but not the transformation 
properties, unless for linear coordinate transformations) are the same at least under some 
conditions. 



VI. GRAVITATIONAL FIELD WITH MATTER SOURCES 

We will now include in our consideration matter fields interacting with the gravitational 
field. One or several matter fields are denoted by where A is some general index. 
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A. Gravitational field equations and the energy-momentum tensor for the matter 

fields 



The total action in presence of the matter sources is 

S= ] -J(L 9 + L m ) d 4 x (79) 



where L 9 is the gravitational Lagrangian ( p4|) and L m is the matter Lagrangian which 
includes interaction of the matter fields with the gravitational field. We assume the universal 
coupling of the gravitational field to all other physical fields, that is, we assume that the L m 
depends on the gravitational field variables h^ v in a specific manner: 



V 



It was shown || that the L m must depend on h^ u and 7 Ml/ only through the combination 
v 7- 7(7^ + h^), if one wants the matter energy-momentum tensor t^ u to participate in the 
gravitational field equations at the equal footing with the gravitational energy-momentum 
tensor, that is, through the total energy-momentum tensor which is the sum of the two. 
The matter energy-momentum tensor t^ v is defined by the previously discussed (see Sec. 
II) universal formula 

2 5L m . . 

^ = ^7 • 81 

V = 7<H^ 

Let us now turn to the derivation of the field equations. The gravitational equations 
are derived by applying the variational principle to the gravitational variables in the total 
Lagrangian. The previously derived equations (^9|) remain unchanged since we assume (for 
simplicity) that the L m does not contain . However, equations ( P5|) are changed and 
take now the form 



r 



Ik 5L m 



0. (82) 



As for the matter field equations, they are derived by applying the variational principle to 
the matter variables in the total Lagrangian, which means: j^- = 0. The concrete form of 
the matter field equations will not be needed. 

We know (Sec. V A) that equations (|8~2"1) without the term caused by L m are equivalent 
to equations ([77]), where Kt^ is given by formula (|65"D. We want to show that the source 
term in the right-hand side of the gravitational equations becomes now, in presence of L m , 
the total energy-momentum tensor. 

Let us start from the contribution provided by L 9 . Since the procedure of reduction of 
Kt^ v to the final form nt^ ( pB] ) involved the use of the equations of motion (|5B"D, which 
are now modified to (|32"D, the gravitational part of the total energy-momentum tensor will 
also be modified, as compared with (|65|) . Using (j82|) instead of (|38D, and getting rid of the 
second derivatives of h^ v in the same way as before, one obtains 



-7 bh a P 
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where the subscript \ m indicates that the derivation has been done in presence of the matter 
fields. The Kt^ u is given of course by the same formula (65), and quantities q a ^^ v are given 



by formula fl59|). Let us now turn to t^ u . The universal coupling in the Lagrangian (|30D, 
that is, the fact that h^ u and 7^ enter the L m only in the combination \/— 7(7^ + h^ u ), 
allows us to relate |^ with J^r. After necessary transformations, one obtains 

1 ST m 



Thus, after using the field equations and removing second derivatives of h^, the total 
energy-momentum 







2 5(L 9 + U 



reduces to 



(7^ + h pv ){^ av + h au ) - -(7^ + h" u ){^ + h af3 ) 



2k 5L m 
/=7 5h a P ' 



Finally, we can write the gravitational field equations in the form similar to equations 
([771). We take the same linear combination of equations (|82|) as was previously done in ( [75|) . 
Putting all the terms in the necessary order, we arrive at the equations equivalent to (|82|): 

I[(y- + hT)^ + h ali ) - (Y a + h^){^ + h vf> )] ;asfi = K(iT\ m + O- ( 83 ) 

Thus, in the gravitational field equations, the total energy-momentum tensor is the source 
for the generalised d'Alembert operator. Obviously, the conservation laws (t fJlU \m + ^ flu );u = 
are satisfied as a consequence of the field equations flB3|) . 

As a final remark, we should mention that the field-theoretical formulation of the general 
relativity allows also gauge transformations in addition to arbitrary coordinate transforma- 
tions. Under gauge transformations, solutions to the field equations transform into new 
solutions of the same equations. In what sense and under which conditions the gauge- 
related solutions are physically equivalent, is a deep and nontrivial issue. This question 
was partially analyzed in ref. |l{J but it is outside of the scope of this paper. The theory 
is fully consistent in its mathematical structure and physical interpretation, if the gauge 
transformations are applied to the gravitational field and matter variables together (even if 
we deal only with a couple of test particles interacting with the gravitational field and which 
are being used in a gedanken experiment). We mention the gauge freedom only in order 
to stress that all the objects and equations have been derived in arbitrary gauge, without 
imposing any gauge conditions. 



B. Equivalence with the geometrical Einstein equations 

In the geometrical approach to the general relativity one interprets the quantities g^, 
introduced by (flop, (f46|), as the metric tensor of a curved space-time (p5|). The universal 
coupling of gravity with matter translates into 
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-ggT\ 



-99 



>A, <PA,a 



(84) 



One can think of this dependence as of manifestation of the Einstein's equivalence principle. 

The matter energy-momentum tensor T^ u is now defined as the variational derivative of 
L rn with respect to what is now the metric tensor: T Mi , = The specific form of the 

L m allows us to write: 



2 5U 



5U 



-99 



pa) 



(85) 



Tensor T^ LU certainly differs from the tensor r^ u defined in the field-theoretical approach, 
but they are related: 



1 



We are now in the position to prove that the field equations ([83]) are fully equivalent to 
the Einstein's geometrical equations. We know (see fl53|) ) that 



R 



R 



(86) 



where R^ = 0. Combining (|82|) and flS5|) we arrive at the Einstein's geometrical equations 



R 



■fat 



« [Tpu - -gpu9 a(S T af3 j . 



(87) 



VII. CONCLUSIONS 



We have shown that the field theoretical formulation of the general relativity allows 
us to derive the fully satisfactory gravitational energy-momentum tensor satisfying all 6 
demands listed in the Abstract of the paper. Both routes, "metrical" and "canonical" , lead to 
one and the same unique expression fl5"5p. When the gravitational field is considered together 
with its matter sources, the same strict rules produce the matter energy-momentum tensor 
t^ v and the modified gravitational energy-momentum tensor. Both tensors participate at the 
equal footing in the nonlinear gravitational field equations ( p3|) which are fully equivalent to 
the Einstein's geometrical equations fl87|) . These strictly defined energy-momentum tensors 
should be useful in practical applications. 
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APPENDIX A: CO VARIANT GENERALISATION OF THE EULER-LAGRANGE 

EQUATIONS 



The covariant field equations (|12"D can be derived in a more traditional fashion, when 
one considers the field variables h^ u and their ordinary (not covariant) derivatives h^ v jT as 
functions subject to variation. To emphasize this fact we rewrite the Lagrangian (|6]) in the 
form 

L = L^\C\ v ^X\ a ) . (Al) 
Starting from the Lagrangian in this form, one derives the usual field equations 

dL ( dL \ 



= . (A2) 

r 

Since the function L in (§) and (|A1[) is one and the same function, but written in terms of 
different arguments, one can relate its derivatives. The second term in ( |A2|) transforms as 
follows: 



dL \ f dL dhP a . ul \ ( dL \ f dL 



\dh^ u T \dhP a . u} dh^ v T \dhP u . T \dh^. T 



+ 



dh av . T ^ dh^, T 
The first term in (|A2|) transforms as 



'"- C'+^C^. (A3) 



dL dL dL dhP a . T dL dL ^ a dL ^ a 



Using ( |A3[ ) and ( ]A4|) in ( [A"2] ) one obtains the required result (|T^). Obviously, the field 



equations ( |A"2| ) use the Lagrangian in the form (|A"T|), whereas the field equations (pT2|) use 
the Lagrangian in the form (M). 



APPENDIX B: PROOF OF EQUATION (68) 



We need to show in detail that the variational derivative of the added term in Lagrangian 
(|66|), calculated at the constraint R ap i3 a = 0, result in the last two terms in equation (|68|). 
Let us introduce a shorter notation for the added term stressing its dependence on 7^ and 
derivatives: 



Taking into account this dependence, we can write for the variation of the function II: 



6U 



dU 



dli 



#7, 



+ 



dli 



flU,T , 



d 7i 
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where 



on 



Slav 



The coefficient in front of 8^ pu defines the variational derivative: 



5U 



<9n 



so that 



fW 



5U 



+ 



dli 



5U 



fJtU 



B 7 ' . 



In order to find the required variational derivative, we can simply calculate the variation of 
■\/—"fA. a P pcr R a pi3 a and present it in the form 

8{y/^A^R apP(T ) = V=^A^6j^ + (V^C T ),r ■ 

The quantity yJ—^A^ is what we need. 

It is convenient to work with yj— r yA af3p ' J 'y aT R r p/3a . instead of y /—'jA al3f " T R ap p (T . The vari- 
ation can be written as 



afipcr , 



.afipo 



7o 



(Bl) 



The first term on the right-hand side of (BT) vanishes due to the constraint, so we need to 
focus attention on the second term. The variation of the Riemann tensor is 



^V) = (^) ; ^-(5^) ;ff . 
The variation of the Christoffel symbols is 



(B2) 



scr^ = i 7 rA (<W + ^A. ; p - s lpa;X ) = \{8Xi Pt + ^i Pt - W-f x )*iaPv ( B3 ) 

One needs to combine (|B~3|), (|B~2| ) and take into account properties of the symmetry of A al3pcr 
(see (|67D). After rearranging the participating terms, one gets 



5(y/^A a ^R apf3a 



[IV 



(B4) 



The second term on the right-hand side of ( p^ ) is a covariant derivative of the vector density, 
so the covariant derivative can be replaced by ordinary derivative, and this term has the form 
of (y/— 7C T ) jT . Thus, we conclude that the sought after variational derivative, calculated at 
the constraint, is 



8{^A^R apPu ) 
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Its contribution to the k t ^ v \ c is — (A M ^ + A Ufial3 ). a .p what we needed to prove. 

The calculation of (|56|) can be done in exactly the same way. Namely, taking the variation 
of the Lagrangian density (^) with respect to 7^, one obtains 



5L 9 



"7 



2k 



2P a pa h pP 8C° afS } . 

Using expression ( p3|) and rearranging the participating terms we arrive at 



(B5) 



5L 9 



-7 



2k 



1 



Thus, as it is stated in the text Q5ED, the variational derivative is 

2k 619 Ktr = lY"h^p%+[>r>r'- 



(B6) 



1 
2 



p P Pa) + g Miy • 



(B7) 
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